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Abstract. New insights into the microscopic origin of itinerant ferromagnetism were 
recently gained from investigations of electronic lattice models within dynamical mean- 
field theory (DMFT). In particular, it is now established that even in the one-band 
Hubbard model metallic ferromagnetism is stable at intermediate values of the inter- 
action U and density n on regular, frustrated lattices. Furthermore, band degeneracy 
along with Hund's rule couplings is very effective in stabilizing metallic ferromag- 
netism in a broad range of electron fillings. DMFT also permits one to investigate 
more complicated correlation models, e.g., the ferromagnetic Kondo lattice model with 
Coulomb interaction, describing electrons in manganites with perovskite structure. 
Here we review recent results obtained with DMFT which help to clarify the origin of 
band-ferromagnetism as a correlation phenomenon. 

1 Introduction 

Since the Curie temperature of ferromagnetic metals like iron, cobalt, and nickel 
is of the order of electrostatic energies in solids, i.e., is much higher than the 
interaction energy of the electron spins, itinerant ferromagnetism is expected to 
be the result of the interplay between the ordinary spin- independent Coulomb 
interaction and the kinetic energy of the electrons, in combination with the Pauli 
principle ■ As such it is one of the fundamental cooperative phenomena in 

condensed matter physics. 

Until recently the theory of itinerant ferromagnetism was investigated by two 
essentially separate communities, one using model Hamiltonians in conjunction 
with many-body techniques (or even rigorous mathematical methods ), 
the other employing density functional theory (DFT) @,||. DFT and its local 
density approximation (LDA) have the advantage of being ab initio approaches 
which do not require empirical parameters as input. Indeed, they are highly suc- 
cessful techniques for the calculation of the electronic structure of real materials. 
However, in practice DFT/LDA is restricted in its ability to describe strongly 
correlated materials. Here, the model Hamiltonian approach is more general 
and powerful since there exist systematic theoretical methods to investigate the 
many-electron problem with increasing accuracy. Nevertheless, the uncertainty 
in the choice of the model parameters and the technical complexity of the cor- 
relation problem itself prevent the model Hamiltonian approach from being a 
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flexible or reliable enough tool for studying real materials. The two approaches 
are therefore complementary. 

Originally the one-band Hubbard model was introduced to gain insight into 



the origin of metallic ferromagnetism |j,|10 11 . However, even for this simplest 
possible microscopic model answers are not easily obtained since in general 
metallic ferromagnetism occurs at intermediate couplings and off half filling 
[p|JT^]. Thus, this cooperative phenomenon belongs to the class of problems 
where standard perturbation techniques are not applicable. In particular, weak- 
coupling theories or renormalization group approaches which are so effective 
in detecting instabilities with respect to antiferromagnetism or superconductiv- 
ity meet with limited success here. In general, non-perturbative investigations 
are required. 

During the last couple of years significant progress was made in the theoret- 
ical understanding of the microscopic foundations of metallic ferromagnetism. 
They were made possible mainly by the development (i) of new analytic ap- 
proaches, such as the mathematical methods used to investigate flat-band fer- 
romagnetism p| and its extensions jl4|,f[5j as well as dynamical mean-field the- 
ory (DMFT) p6p7| , p^[ |, (ii) of new numerical techniques, such as the density 
matrix renormalization group (DMRG) which yields precise results in d = 1 
[E9I, and (iii) of new comprehensive approximation schemes such as the multi- 
band Gutzwiller wave function or the new ab initio computational scheme 
LDA+DMFT |2l|,||,|3) which merges conventional band structure theory (LDA) 
with a comprehensive many-body technique (DMFT). 

In this paper the insights gained with the help of the DMFT will be reviewed. 
After a short introduction into this approach (Section , the microscopic con- 
ditions for metallic ferromagnetism in the one-band Hubbard model (Section ||) 
and in the case of the orbitally degenerate model (Section |j) are explained and 
the differences discussed. Furthermore, the physics of itinerant ferromagnetism 
in more complicated models, e.g., the ferromagnetic Kondo lattice model with 
Coulomb correlations for manganites with perovskite structure, is analyzed. A 
conclusion and outlook (Section ||) closes the presentation. 



2 Dynamical Mean-Field Theory 

During the last decade a new many-body approach was developed which is es- 
pecially well-suited for the investigation of correlated electronic systems with 
strong local interactions - the dynamical mean-field theory (DMFT) [Q. It 
becomes exact for d — 00, i.e., for lattices with coordination number Z = 00 
[ fi6| . Why is the investigation of this limit useful and interesting? The answer 
is that already in d = 3 the coordination number of regular lattices, such as 
the fee lattice, is quite large (Zf cc = 12). It is therefore quite natural to view Z 
as a large number, and to consider the limit Z — > 00 as a starting point fjq| . 
To obtain a meaningful model in this limit one has to scale the NN hopping 
amplitude in the kinetic energy (see below) as t = t* /yZ (in the following we 
set t* — 1 and thereby fix the energy scale). Then one obtains a purely local 
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theory where the self-energy Sk{u>) becomes k independent and where the prop- 
agator G(k, u>) — G°(k, u> — 2J{ui)) may be represented by the non-interacting 
propagator at a shifted frequency |24| . The dynamics of the quantum mechanical 
correlation problem is then fully included. That is why this theory is referred 
to as "dynamical mean- field theory" (DMFT). We note that due to the local 
nature of the theory there is no short-range order in position space. 

Within DMFT the electronic lattice problem is mapped onto an effective self- 
consistent single-impurity Anderson model |2^| . DMFT is a non-perturbative, 
thermodynamically consistent theoretical framework within which the dynam- 
ics of correlated lattice electrons with local interaction can be investigated at 
all coupling strengths. This is of essential importance for problems like band- 
ferromagnetism or the metal-insulator transition which set in at a Coulomb 
interaction strength comparable to the electronic band width. 

In DMFT the information about the lattice or the dispersion of the system 
under investigation enters only through the density of states (DOS) N°(E) of 
the non-interacting particles, unless there is long-range order with broken trans- 
lational symmetry of the lattice as in the case of antiferromagnetism. In finite 
dimensions, e.g., d = 3, DMFT is then an approximation - usually an excellent 
one as is manifested by the plethora of results obtained within the last decade 
[ p4| . DMFT is presently the only comprehensive, thermodynamically consistent 
computational scheme which allows one to investigate the dynamics of correlated 
lattice electrons on all energy scales. 

Due to its equivalence to an Anderson impurity problem a variety of ap- 
proximative techniques have been employed to solve the complicated DMFT 
equations, such as the iterated perturbation theory (IPT) |2(| and the non- 
crossing approximation (NCA) |27|, as well as numerically exact techniques like 
quantum Monte-Carlo simulations (QMC) p8| , exact diagonalization (ED) plfl , 
or the numerical renormalization group (NRG) |50|]. However, NRG cannot yet 
be used to solve the DMFT equations for multi-band models. 

In the present paper all DMFT results were obtained by QMC. This is a 
particularly well-tested and reliable, albeit computer-expensive method which 
may be employed down to temperatures T ~ lO -2 !^ (W: band width) and at 
not too strong interactions. 



3 The One-Band Hubbard Model 

The one-band Hubbard model 

#Hub = - ( C L C 3* + h - C -) + U Y1 n it n U C 1 ) 

is the simplest lattice model for correlated electrons, and was originally proposed 



to understand metallic ferromagnetism in 3d transition metals |9|,[10 11 . Here, t%j 
is a general hopping matrix element between sites i and j. In the past the kinetic 
energy in ([!]) was usually restricted to nearest-neighbor (NN) hopping; then it 
is useful to divide the underlying lattices into bipartite and non-bipartite ones. 
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Fig. 1. DOS of noninteracting electrons on a fee lattice in d = 3 with and without 
additional NNN hopping t' . 



This distinction ceases to be useful if, in addition to NN hopping t, longer-range 
hopping, in particular next-nearest-neighbor (NNN) hopping t' is considered. 

The Hubbard model is characterized by a purely local interaction term which 
is completely independent of the lattice structure and dimension. Therefore in 
this model the kinetic energy of the electrons (i.e., the dispersion) and the struc- 
ture of the lattice are very important for the stability of metallic ferromagnetism. 
This is well-known from the approximate investigations by Gutzwiller Q , Hub- 
bard |l0|, Kanamori pf and Nagaoka pM, and has recently been confirmed, 



and made precise, by detailed investigations P,p2Ula,14Ll / 18 



3.1 Routes to Ferromagnetism 

On bipartite lattices the i'-hopping term destroys the antiferromagnetic nesting 
instability at small U (see, for example |33|) by shifting spectral weight to the 
band edges and thereby introducing an asymmetry into the otherwise symmetric 
DOS. It will be shown below that a high spectral weight at the band edge (more 
precisely: the lower band edge for n < 1) minimizes the loss of kinetic energy of 
the overturned spins in the magnetic state and is hence energetically favorable. 
Therefore frustrated, i.e., non-bipartite lattices, or bipartite lattices with frus- 
tration due to further-range hopping (e.g., if ^ 0) support the stabilization of 
metallic ferromagnetism. 

The fee lattice is an example of a frustrated lattice in d = 3. The corre- 
sponding DOS of the non-interacting particles is shown in Fig. [|. Switching on 
an additional NNN hopping t' further increases the spectral weight at the lower 
band edge. For t' = t/2 one even obtains a square-root-like divergence. To un- 
derstand why a high spectral weight at the band edge is favorable for the kinetic 
energy we first consider the case U = 0, n < 1 Jl^]. Let us consider a flat, sym- 
metric DOS. The fully polarized state is obtained by inverting the spin of the 
down electrons, which due to the Pauli principle have to occupy higher energy 
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states. Counting the energy from the lower band edge the Fermi energy of the 
polarized state, fj,^, is seen to be twice that of the unpolarized state. This should 
be contrasted with a DOS having large spectral weight at the lower band edge. 
Here the Fermi level of the polarized state is not so strongly shifted upwards, 
i.e., fewer high energy states are populated, which is clearly energetically favor- 
able. The energy difference between the fully polarized state and the unpolarized 
state 
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must become negative for the ferromagnetic state to be stable. Of course, in the 
non-interacting case one has AE > Q. Nevertheless, even for U = 0, AE 
attains its lowest value for a DOS with peaked spectral weight at the lower band 
edge for all n < 1 [^2|. To show that AE < for U > requires a reliable 
calculation of the energy of the correlated paramagnet. Indeed, this is one of 
the central, most difficult problems of electronic correlation theory. It should be 
noted that the above discussion concerning the shape of the DOS goes beyond 
the well-known Stoner criterion which predicts a ferromagnetic instability of the 
paramagnet for U equal to the inverse of the DOS precisely at the Fermi level. 

3.2 Frustrated Lattices 

Since metallic ferromagnetism is an intermediate-coupling problem purely an- 
alytic approaches alone cannot provide sufficient information, in particular in 
dimensions d > 1. In this situation the development of new numerical tech- 
niques during the last few years was of crucial importance for progress in this 
field. In particular, for one-dimensional systems the DMRG p9[ , and for high- 
dimensional systems the DMFT have led to important insights. Here we restrict 
ourselves to results obtained by DMFT. 

By a suitable generalization of the dispersion relation in three dimensions, 
frustrated lattices like the fee lattice can be defined in any dimension, in partic- 
ular in d = oo [j34| . With the proper scaling of the hopping term (see above) the 
non-interacting DOS of the generalized fee lattice in d — oo is obtained as: 

N° gfcc (E) = e-^ E V*/y/*{l + y/2E) (3) 

It has a strong square- root divergence at the lower band edge, — 1/ V2, and no 
upper band edge. 

Investigations of the stability of metallic ferromagnetism on fec-type lattices 
within DMFT, in combination with finite-temperature QMC techniques to solve 
the DMFT equations, were first performed by Ulmke |p5| . To detect a ferromag- 
netic instability one calculates the temperature dependence of the uniform static 
susceptibility \f from the corresponding two-particle correlation function |3(| . 
For Ngj cc (E) and at an intermediate interaction strength of U = 4 the ferromag- 
netic response is found to be strongest around quarter filling (n ~ 0.5). The sus- 
ceptibility xf is seen to obey a Curie- Weiss law (Fig. ^]). Thus the Curie temper- 
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Fig. 2. Magnetization m (diamonds) and inverse ferromagnetic susceptibility Xf* 
(circles; note the factor of 4 in the definition) as calculated by DMFT for the one-band 
Hubbard model with DOS (|^) at U = 4 and n = 0.58. Error-bars are of the size of the 
symbols or smaller. The value of Xf 1 at T = 0.05 is a data point, not an extrapolation. 
The dotted line is a linear fit to Xf 1 ' ^ ne dashed line a fit with a Brillouin function to 
m H§. 



ature T c can be safely extrapolated from the zero of Xf ■ For n = 0.58 this leads 
to a value of T c — 0.051(2). Below T c the magnetization m (a one-particle quan- 
tity determined by the local propagator) grows rapidly, reaching more than 80% 
of the fully polarized value (m max = n = 0.58) at the lowest temperature (30% 
below T c ). The three data points m(T) in Fig. |^ are consistent with a Brillouin 
function with the same critical temperature of T c = 0.05 and an extrapolated full 
polarization at T = 0. (We note that a saturated ground-state magnetization 
is also consistent with the single-spin-flip analysis of the fully polarized state 
by Uhrig 17 ) . So one finds simultaneously a Curie- Weiss-type static suscepti- 
bility with Brillouin-function-type magnetization, and a non-integer magneton 
number, in qualitative agreement with experiments on 3d transition metals. In 
the past, these features were attributed to seemingly opposite physical pictures: 
the former to localized spins, and the latter to itinerant electrons. However, this 
was only because of the use of oversimplified, static mean-field-type theories 
like Weiss mean- field theory (for spin models) and Hartree-Fock (for electrons). 
Here we see that these properties are natural features of correlated electronic 
systems, which are generated by the quantum dynamics of the many-body prob- 
lem. Within DMFT this "paradoxical" behavior of band ferromagnets is resolved 
without difficulty. 

Collecting the values of T c (U,n) obtained by x^i^c U,n) = 0, the bound- 
aries between the ferromagnetic and the paramagnetic phases are determined. 
Thereby one can construct the T vs. n phase diagram for different values of U. 
The region of stability increases with U j35| . 

To make contact with d — 3 one can use the corresponding fee DOS shown 
in Fig. |l|. For t' = no instability is found at temperatures accessible to QMC. 
However, already a small contribution of <'-hopping (which is present in any 
real system) is enough to stabilize a large region of metallic ferromagnetism in 
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Fig. 3. (a) T vs. n phase diagram of the Hubbard model as obtained within DMFT 
for the DOS corresponding to a three-dimensional fee lattice with NN hopping t' = i/4 
(see Fig. Q); (b) DOS of the interacting electrons in the ferromagnetic phase of (a), 
solid line: majority spin, dashed line: minority spin |3q]. 



addition to an antifcrromagnetic phase (which is absent in d = oo), close to 
half filling (Fig. [|a) |55) . This shows the strong and subtle dependence of the 
stability on the dispersion and the distribution of spectral weight in the DOS. 
By developing increasingly refined schemes of iterated perturbation theory to 
solve the DMFT equations Nolting and collaborator |37| obtained T vs. n phase 
diagrams which reproduce the QMC results remarkably well. 

We note that the maximal transition temperature is T ( ! nax = 0.05 ~ 500 K 
for a band width W = 4 eV. This is well within the range of real transition 
temperatures, e.g., in nickel. 

So far we only argued on the basis of the shape of the DOS of the non-in- 
teracting electrons, N°(E). On the other hand the interaction will renormalize 
the band and relocate spectral weight. Therefore it is not a priori clear at all 
whether the arguments concerning the kinetic energy discussed in subsection 3.1 
still hold at finite U. We note that this effect is also not taken into account in 
any Stoner theory. To settle this point one may calculate the DOS of the inter- 
acting system, N(E), by employing the maximum entropy method for analytic 
continuation. In Fig. ||b we show N(E) obtained with N^ CC {E) and the param- 
eter values used to calculate Fig. ||. Clearly the ferromagnetic system is metallic 
since there is appreciable weight at the Fermi level (E = /i). Furthermore, the 
spectrum of the majority spins is seen to be affected only slightly by the interac- 
tion, the overall shape of the non-interacting DOS being almost unchanged (the 
magnetization is quite large (m — 0.56 at n = 0.66) and hence the electrons in 
the majority band are almost non-interacting). This implies that the arguments 
concerning the distribution of spectral weight in the non-interacting case and 
the corresponding kinetic energy are even applicable to the polarized, interact- 
ing case. The spectrum of the minority spins is slightly shifted to higher energies 
and has a pronounced peak around E — /i ~ U = 6. 

From Fig. |^b we find the exchange splitting A between the majority and 
minority bands to be A = 0.8. Comparing the quantity Um {— 3.4 in the 
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Fig. 4. (a) T vs. n phase diagram of the Hubbard model as obtained within DMFT; 
(b) corresponding shapes of the non-interacting DOS; Fermi energies for n = 0.3 are 
indicated by vertical lines [B8|. 



present case) with A and T c we find a characteristic hierarchy of energy scales: 

Urn > A > T c . (4) 

This is very different from results obtained by Stoner theory or effective one- 
particle theories like LDA where all three quantities are essentially equal. The 
generation of small energy scales is a genuine correlation effect. 

To study the influence of the distribution of spectral weight on the stability 
of ferromagnetism within the DMFT systematically, Wahle et al. [Q recently 
solved the DMFT equations with a tunable model DOS, 



with c = (l+y/T—ar) / (jrD) and half-band width D = 2. Here a is an asymmetry 
parameter which can be used to change the DOS continuously from a symmetric, 
Bcthc lattice DOS (a = 0) to a DOS with a square-root divergence at the 
lower band edge (a = 1), corresponding to a fee lattice with t' = t/A in d = 3 
(Fig. [I]). The strong dependence of the stability of metallic ferromagnetism on 
the distribution of spectral weight is shown in Fig. ^a. Already a minute increase 
in spectral weight near the band edge of the non-interacting DOS, obtained by 
changing a from 0.97 to 0.98 (see Fig. ||b) is enough to almost double the stability 
region of the ferromagnetic phase. Obermeier et al. |39|| found ferromagnetism 
even on a hypercubic, i.e., bipartite, lattice, but only at very large U values 
(U > 30). 

The importance of genuine correlations for the stability of ferromagnetism is 
apparent from Fig. |5[ where the DMFT results are compared with Hartree-Fock 
theory pSfl . The quantum fluctuations, absent in Hartree-Fock theory, are seen 
to reduce the stability regime of ferromagnetism drastically. Spatial fluctuations 
(e.g., spin waves), absent in the DMFT, may reduce the stability regime further. 
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Fig. 5. T vs. U phase diagram for a strongly peaked DOS (a = 0.98, see Fig. ^b) at 
n — 0.4 within DMFT (circles; dashed line is guide to the eyes only) in comparison 
with Hartree-Fock (solid line) [Q. 



3.3 Additional Interactions 

In the one-band Hubbard model only the local interaction is retained. Thereby 
several interactions which naturally arise when the Coulomb interaction is ex- 
pressed in Wannier representation are neglected. Even in the limit of a single 
band and taking into account only NN contributions, four additional interactions 
appear pp|Mpl: 



(0) 



NN L a 

-2F(Si ■ Sj + jTHrij) + F'(4i c li c ji c j] + h - c -) 



Here the first term corresponds to a density-density interaction, the second term 
to a density-dependent hopping, and the fourth term describes the hopping of 
doubly occupied sites. In particular, the third term (with F — F* /Z > 0) 

ff ir = _2^^S i -S j (7) 

NN 

describes the direct ferromagnetic exchange between electrons on NN sites. It is 
this interaction which Hciscnbcrg singled out in his original spin-model as the 
main source of ferromagnetism. It should be noted, however, that this interaction 
is present even when the electrons are free to move. The exchange interaction 
will be quite small, but nevertheless it favors ferromagnetic ordering in the most 
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Fig. 6. Bond-charge interaction X vs. U phase diagram for the generalized model 
H — -Hfjub + Hx within the Gutzwiller approximation for the semi-elliptical Bethe 
DOS at n = 0.9 @. 



obvious way. Hirsch J4l| argued that this term is the main driving force for 
metallic ferromagnetism in systems like iron, cobalt, and nickel. Indeed, one can 
show rigorously that a next-neighbor direct exchange interaction, if only chosen 
large enough, easily triggers the ferromagnetic instability |4^ , |l^ , ^8| . Hence the 
NN exchange may well be important for systems on the verge of a ferromagnetic 
instability. 

Another term that can be of importance for the stabilization of ferromag- 
netism is the bond-charge interaction (X = X* j\[Z) 



Hx = 



X 



^ NNct 



hx.)^ 



+ Uj 



0- 



(8) 



It effectively gives rise to correlated hopping, i.e., the hopping amplitude of an 
electron now depends on the presence of electrons of the opposite spin polariza- 
tion. The magnitude of X has been estimated to be of order 0.1-1 eV pC|pO[ |, 
and hence X may be comparable to the hopping t, although typically smaller 
than U. Hartree-Fock treatment of H — -ffnub + Hx shows that correlated 
hopping can lead to a spin-dependent band narrowing which may stabilize ferro- 
magnetism 1 43 . We note that in DMFT Hx does not reduce to its Hartree-Fock 
contribution. The correlation effects introduced by Hx were recently studied 
within the Gutzwiller approximation (GA) ]l4|. The GA yields the exact evalu- 
ation of expectation values in terms of the Gutzwiller wave function in the limit 
of d — > oo [ 45|Jl^1 and goes beyond Hartree-Fock theory by including correlations 
explicitly. It was found that correlated hopping with X > can lower the crit- 
ical value of U for ferromagnetism considerably as shown in Fig. ^. Compared 
to Hartree-Fock theory the ferromagnetic region of the phase diagram is seen 
to be much reduced. Recently it was also shown how to incorporate Hx into 
DMFT at least in principle. It will be interesting to see how the GA results 
compare with future DMFT calculations for H = i/nub + Hx- 



Metallic Ferromagnetism - an Electronic Correlation Phenomenon 11 



Clearly, NN interactions may help to stabilize ferromagnetism. However, since 
the considerably larger Hubbard interaction U, together with a suitable kinetic 
energy, is already sufficient to trigger a ferromagnetic instability, the ferromag- 
netic exchange and bond-charge interactions appear to play only a secondary 
role. 



4 Band Degeneracy and Local Exchange 

Besides the NN Heisenberg exchange interaction another much larger exchange 
term is present in ferromagnets like iron, cobalt, and nickel, namely the local 
exchange between electrons in different orbitals on the same lattice site. It has 
long been speculated that this exchange interaction, which is known to align 
electrons on isolated atoms (Hund's first rule), may also lead to ferromagnetism 
in the bulk, being transmitted by the kinetic energy J47| . The simplest model 
for this mechanism is the two-band Hubbard model with local exchange Fq and 
Coulomb repulsion Vq between two orbitals v = 1, 2 (Fig. |7|): 

#2-band = -* c tv<r c jv<r + u 5Z n M"M 

NN, av iv 

+ J2(V - (W'^oKio-rcaw - p aJ2 c ii* c iia> c i2a> c i2a- ( 9 ) 

iaa' a^a' 

This two-band model Kj| and its multi-band generalizations H] were studied in 
considerable detail by various theoretical techniques p(J. Recent investigations 
[ pl| were triggered by the renewed interest in the electronic properties of transi- 
tion metal oxides where the doubly degenerate e g bands of the d electrons plays 
a very important role. At quarter filling (one electron per site; n — 1), ferromag- 
netism can be understood by superexchange within strong-coupling perturbation 
theory (Fig. |J) . But away from quarter (or half) filling, the model is much more 
difficult to treat. In this regime DMFT, solved by QMC, once more provides a 
powerful tool to investigate this model |p2fl , at least, if the last term in @, i.e., 
the spin-flip contribution of Fq , is neglected. 

The calculated T vs. n phase diagram is shown in Fig. [| for a symmetric 
Bcthc DOS. Without Hund's rule coupling Fq ferromagnetism is unstable for 



this DOS, at least at moderate values of U (see Sec. |3.2| for the one-band model 
and (5^] for the two-band model). However, when the Hund's rule exchange 
is included a metallic ferromagnetic phase is stabilized between an alternating 
(staggered) orbital ordered state at quarter filling and antiferromagnetism at 
half filling. 

While orbital ordering at n = 1 and antiferromagnetism at n = 2 can be 
understood by superexchange, the mechanism for metallic ferromagnetism is 
more subtle. Of course, the virtual superexchange processes described above 
are still present. But the additional electrons doped to the n = 1 system may 
now move freely between two singly occupied sites, i.e., are unhindered by the 
Coulomb repulsion U or Vq. Furthermore, if the electrons on the two sites are 
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Fig. 7. Illustration of the local interactions (Hub- 
bard U and Hund's rule couplings Vo, Fo defined in 
(0)) between electrons in a two-band model. 



4 4 

~2t 2 /U -2t 2 /V -2t 2 /(V -F ) 

Fig. 8. Energy gain in second-order perturbation theory in t for four configurations 
with two electrons on two sites. Most favorable is ferromagnetism combined with 
an alternating occupation of orbitals. Note, that the true ground state is not the 
last configuration above, i.e., ct_ lj/=:1 ^c^ a « |0) , but the corresponding orbital singlet 

2( C 11| C 22T — C 12 ^C 21 ^)|0). 



AE= 



spin aligned they do not even have to spend the exchange energy Fq. Therefore a 
ferromagnetic alignment of the spins improves the kinetic energy of the electrons. 
This is essentially the concept of double exchange introduced by Zener p3| to 
explain ferromagnetism in manganites such as Lai-^CaxMnOs, and put on a 
firmer theoretical fundament by Anderson and Hasegawa [ pl[ . 

In manganites, the cubic crystal field splits the five Mn d orbitals into three 
t2 g and two e g orbitals. The former have a lower energy and hybridize less 
strongly with the O p orbitals. Thus, the three electrons within the t^g or- 
bitals can be approximately described by a localized spin Si , with the remaining 
1 — x electrons occupying the e g orbitals. If the Coulomb repulsion between e g 
electrons is taken into account one arrives at a correlated electron model for 
manganites: 

2 

H = i?2-band — 2 J ^ S iL , ■ S^. (10) 
v— 1 i 

Here, i/ 2 -band is defined in (|), s iv =\ Xw c \v<7 T <?<?' c iv<j> denotes the e g spin (r: 
Pauli matrices), and J is the local exchange between t2 g and e g electrons. With- 
out Coulomb interaction (U — Vo = Fq = 0), the Hamiltonian ( |Io| ) reduces 
to the ferromagnetic Kondo lattice model (KLM) which was investigated inten- 
sively in recent years ]5q] . This model forms the theoretical basis for the double 
exchange mechanism: at J > f, the optimization of the kinetic energy of the 
e g electrons requires a ferromagnetic environment of t2 g spins. The KLM fails 
to describe the correct behavior for x < 0.5 since it does not penalize double 
occupations, i.e., two e g electrons on the same site. The suppression of double oc- 



Metallic Ferromagnetism - an Electronic Correlation Phenomenon 13 




Fig. 9. T vs. n phase diagram of the two-band Hubbard model including antiferro- 
magnetism (AF), ferromagnetism (FM), and orbital ordering (OO) for a Bethe DOS 
(total width W = 4), U = 6, Vo = 4, and F = 2 §§. 



0.15 



T 0.1 




0.05 



Fig. 10. Curie temperature T c for the phase transition from the paramagnetic (PM) 
to the ferromagnetic (FM) phase as a function of Ca doping x (n = 1 — %). Dashed 
line: KLM with Bethe DOS (band width W = 2) and Hund's rule coupling J = 3/2; 
squares: correlated electron model ( |lo[ ) which also takes into account the Coulomb 
interaction between e g electrons (U = 8, Vo = 6, and Fq = 1); cross: Weiss mean-field 
theory (superexchange) for this model [M. Note that, without the coupling to the t^g 
spin, no ferromagnetism is observed in Fig. ^ for n < 1 at similar interaction strengths. 



cupations induced by the Coulomb repulsions U and Vo was investigated within 
DMFT in [f36f , and was shown to lead to a crossover from double exchange to 
superexchange. This results in a maximum in the Curie temperature (Fig. |Lo| ) 
in qualitative agreement with experiment |57j . 
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5 Conclusion and Outlook 

In this paper we reviewed recent developments in our understanding of the ori- 
gin of metallic ferromagnetism in the one-band Hubbard model and in band- 
degenerate models on the basis of the dynamical mean- field theory (DMFT). 
In the one-band Hubbard model metallic ferromagnetism is found to exist in a 
surprisingly large region of the U vs. n phase diagram. A stabilization of this 
phase at intermediate U values requires a sufficiently large spectral weight near 
one of the band edges as is typical for frustrated (e.g., fee- type) lattices which 
optimize the kinetic energy of the polarized state and, at the same time, frustrate 
the parasitic antiferromagnetic order. This finding, together with the results ob- 
tained for dimension d = 1 |^ , [i^ , |l4|Jl9| , finally establishes the stability of band 
ferromagnetism in the one-band Hubbard model on regular lattices and at inter- 
mediate values of the interaction U and density n. Thereby one of the prominent 
questions of many-body theory in this field is finally answered. 

By contrast, the origin of metallic ferromagnetism in the band-degenerate 
Hubbard model at intermediate U values is not primarily a spectral weight effect 
but is already caused by moderately strong Hund's rule couplings . In this respect 
the origin of ferromagnetism in the orbitally degenerate model is more straight- 
forward than that in the (less realistic) one-band case. Nevertheless, in the ab- 
sence of orbital ordering the resulting magnetic phase diagrams are remarkably 
similar. The identification of a single main driving force for the stabilization 
of metallic ferromagnetism in the one-band and the band-degenerate Hubbard 
model, respectively, helps to differentiate between different mechanisms. How- 
ever, in real systems these effects will tend to conspire, as is evident, for example, 
in nickel where an fee lattice leads to a strongly asymmetric DOS and the band 
degeneracy brings with it Hund's rule couplings. 

These insights also helped to understand the origin of itinerant ferromag- 
netism in more complicated models, e.g., the ferromagnetic Kondo lattice model 
with and without Coulomb correlations, which is employed to understand the 
physics of manganites with perovskite structure, like Lai_ a; (Sr,Ca) 2 ;Mn03. It 
was found that in this model double exchange can explain ferromagnetism only 
for doping x > 0.5. At lower values of x the suppression of double occupations 
by the local Coulomb repulsion becomes more and more important and leads to 
a crossover from double exchange to superexchange. This results in a maximum 
of the Curie temperature in qualitative agreement with experiment. 

As discussed in the Introduction the problem of metallic ferromagnetism was, 
until recently, investigated by two essentially separate communities - one using 
model Hamiltonians together with many-body techniques, the other employing 
density functional theory (DFT) in the local density approximation (LDA). In 
view of the individual power of LDA and the model Hamiltonian approach it 
is highly desirable to combine these techniques, thereby creating an enormous 
potential for all future ab initio investigations of real materials, including, e.g., 
/-electron systems, Mott insulators and metallic ferromagnets. A combination 
of these two approaches had already been used to investigate band ferromag- 
netism some time ago fl58| . Recently, a fusion has started to emerge in new 
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directions. One is the construction of multi-band Gutzwiller wave functions in 
combination with spin-density functional theory in the limit of large coordina- 
tion numbers p0| . It was already successfully applied to ferromagnetic nickel, 
leading to significant improvements over LDA results. The other is the ab initio 
computational scheme LDA+DMFT |23] which was recently used to investigate 
transition metal oxides and to calculate the magnetic excitation spectrum of fer- 
romagnetic iron . Without doubt these and related methods [^9j will rapidly 
develop into standard tools for future investigations of band ferromagnetism and 
other electronic correlation phenomena. 
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